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A model of carrier phase measurement (as carried out by a satellite navigation receiver) is 
04 ' formulated based on electromagnetic theory. The model shows that the phase of the open-circuit 

\ voltage induced in the receiver antenna with respect to a local oscillator (in the receiver) depends 

, on the relative orientation of the receiving and transmitting antennas. The model shows that using 

a single receiving antenna, and making carrier phase measurements to seven satellites, the 3-axis 
O , attitude of a user platform (in addition to its position and time) can be computed relative to an initial 

point. This measurement model can also be used to create high-fidelity satellite signal simulators 
that take into account the effect of platform rotation as well as translation. 
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I. INTRODUCTION 



' O , Satellite systems such as the Global Positioning System (GPS) and the Global Navigation Satellite System 
^ ' (GLONASS) were originally designed to provide a user with position and time . In addition to position and time, 

^ , orientation (with respect to some system of coordinates) or attitude information has numerous applications, such as 
to aircraft flying by instrument navigation, pointing communications antennas and satellites, unmanned air vehicles 
(UAV), and as a source of compass heading. In recent years, attitude (or orientation) information has been obtained 
from the GPS by using multiple, widely-separated antennas mounted on a (assumed) rigid platform When the 
rigid platform translates and rotates, the path length changes between the transmitting satellite and each receiving 
(-^ antenna. Each satellite transmits a (modulated) continuous wave sinusoidal signal. The phase of the voltage induced 
O i' in each receiving antenna changes with respect to the phase of a local oscillator in the receiver. A measurement of 
this effect is called a carrier phase measurement (or delta-range measurement) and is usually modeled Q as the phase 
difference (or accumulated phase difference from some point in time) between the receiver's internal oscillator and the 
^ \ voltage induced on the given receiving antenna. 

QQ I In current GPS receiver technology, this phase change is modeled as arising solely from the change in path length 
04 ' between transmitting satellite and receiver antenna [Q. This model is used in commercial GPS receivers to determine 
. the current orientation or attitude of the rigid platform (on which the antennae are mounted) relative to the (platform) 
attitude at some initial point in time. The standard model of carrier phase measurement Q completely neglects the 
physical nature of the electromagnetic field (i.e., that it is a vector field) and the nature of the detection process that 
occurs in a real antenna. The role of the antenna is to convert the vector electromagnetic field into a time-dependent 
voltage signal. 

^ In recent work, Krall and Bahder [^,^ developed an improved physical model of the carrier phase measurement 

• I— I [ process. They have shown that using a single dipole antenna two angles of attitude can be obtained by tracking six 
satellites (|,^. Their model of carrier phase measurement treats the electromagnetic field as a vector and includes 
(— I ' orientation effects of the detection process of the receiving antenna. In particular, the model of Krall and Bahder 
^ y shows that the phase of the open circuit voltage induced in the receiving antenna depends on the relative orientation 
!| ■ of the receiving and transmitting antennas, which is a feature that is neglected by the standard model Q of carrier 
. ^ ' phase measurement . 

The equations presented by Krall and Bahder ||,^ are specific to a dipole receiving antenna and these equations 
?H cannot be simply extended to more realistic types of receiving antennas, such as those that are currently used in actual 
GPS receivers. Furthermore, the Krall and Bahder equations allow determination of only two angles of attitude, while 
many applications require knowledge of all three angles of attitude for a user platform. 

In this paper, I generalize the work of Krall and Bahder. I show that, for a general type of antenna, all three angles 
of attitude can be computed from carrier phase measurements, i.e., when the receiving antenna has no cylindrical 
symmetry. I present a generalized model of the carrier phase measurement process based on the vector effective length 
of an arbitrary antenna. The vector effective length characterizes a receiving antenna in terms of a two-component 
complex vector field, h(6', </>), on the surface of a unit sphere surrounding the antenna [^,| 10 1^. The spherical polar 



angles 0, and (j) specify the direction from which the radiation is incident on the receiving antenna, as seen in the 
comoving frame of the receiving antenna. 
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The improved model of carrier phase measurement developed below can be used to compute all three angles of 
platform attitude from single- antenna carrier phase measurements to seven satellites. Since the phase of the voltage 
induced in the receiving antenna depends on the relative orientation of the transmitting and receiving antennas, the 
attitude of the satellite antennas must be known to the user of the system. For the case of GPS, the satellite orientation 
is quite strict because the power in the main antenna beam is centered on the Earth [T2| . So the satellite antenna 
attitude is assumed known in the model presented below. However, for future high-accuracy attitude determination, 
the satellite system should have its navigation message enhanced so that it contains information on the orientation of 
each satellite antenna. As a result of this enhancement, the GPS would be a a navigation system that provides users 
with their position, time, and orientation (or attitude). 

The single- antenna feature of the model presented below permits miniaturization of the attitude determination 
system over that of present systems, since multiple antennas are not needed. Consequently, this model is advantageous 
in applications that require miniaturization, such as pointing communications antennas, nano-satellites JTsf , micro air 
vehicles (MAV) and wrist-mounted compass heading systems This model can also be used for developing 
a satellite signal simulator that faithfully models carrier phase measurements made by a platform that undergoes 
rotation (in addition to translation) such as a satellite, aircraft or missile. 

Additionally, in the currently-used multiple-antenna attitude determination systems, higher accuracy of attitude 
is achieved with larger antenna separations (baselines). However, in systems with large baselines, alignment of 
the antennas, as well as change of antenna orientation due to flexing of the (assumed rigid) mounting base, are 
complications that can lead to significant phase errors. For this reason, a single-antenna attitude determination 
system may be preferable. Alternatively, it may be advantageous to use the model of carrier phase measurement 
presented below in multiple antenna systems so that the requirements of platform rigidity and antenna alignment can 
be relaxed. 

The apparent phase center of a receiving antenna moves as a function of direction of the incident radiation from 
the satellite. The direction from which the signal arrives at the receiving antenna, as seen in the comoving frame 
of the receiving antenna, depends on the relative velocity of the transmitting and receiving antennas (aberration of 
star light effect). This effect was taken into account by Krall and Bahder to first order in v/c '--^lO"^, where v is the 
speed of the satellite antenna in the comoving frame of the receiving antenna, and c is the speed of light. For typical 
satellite navigation applications, the speed of the receiver is small (compared to c) and this aberration effect amounts 
to approximately 10~^ radians of angle, and therefore can be neglected. In this work, I neglect this first order effect 
by dropping terms of first order in v/c. However, since the orientation effects (such as dependence of the phase of the 
voltage on direction of incident radiation and movement of the apparent phase center) are zeroth-order effects, they 
are included in the model below. 

In this paper, the standard model of carrier phase measurement is reviewed in section II. Section HI, contains a 
calculation of the open-circuit voltage induced in a receiving antenna due to the satellite signal. The carrier phase 
measurement model based on the open-circuit antenna voltage is formulated in section IV. Section V deals with 
application of the measurement model to attitude determination and navigation. Section IV presents a summary of 
the ideas involved in this work. 



The standard multiple-antenna technique of attitude determination relies on a simplified model of carrier phase 
measurements [Q. This standard technique makes use of differences between the phase of a local oscillator and the 
phase of the incoming satellite signal, at the positions of the (multiple) receiving antennas. In this technique, the 
satellite signal is essentially treated as a scalar wave of amplitude E, given by 



where Eo is a real amplitude, r = |r — R| is the distance from the satellite antenna at transmission point R to receiver 
antenna at reception point r, the angular frequency co = 27r/, / is the frequency of the carrier signal and j = \/"-T. 
The wave propagates with wave vector k = Nui/c — 2ttN/X, at speed c/N, where c is the speed of light and N is 
the average index of refraction over the signal's path. In the simplified model of carrier phase measurements, the 
measured phase is taken to be the phase of the scalar field, seen by antenna i at time t at position ri(t), due to a 
satellite j at position Rj {t — t): 



II. STANDARD MODEL OF CARRIER PHASE MEASUREMENT 



E = Eo e-J('=''-"*) 



(1) 



(j)^j{t) ^ k\r,{t) ~ Rj{t - t)\ - ujt 



(2) 



where r is the time of flight of the signal, implicitly defined by 
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^ ^ |r,ft)-Rj-(t-r)| 

c 

The satellite receiver has a local oscillator that is not perfect. The imperfection of the oscillator is represented by 
the fact that the oscillator does not keep time according to Earth-centered inertial (ECI) coordinate time or universal 
coordinated time (UTC), t, but according to its own clock that counts cycles according to a time scale t* , which 
is some function of, say, ECI time t. The oscillator will then have a running phase given by 

(boUt) = ~^t* (4) 

where the time scale t* is not the same as the ECI (or UTC) time scale t, but there is a functional relation t* = t*{t) 
between these two times. The change in phase of the local oscillator since some initial time or epoch, to, can then be 
represented by 

Acjjoscit) = -uj{t* - to) (5) 

The difference between the ECI coordinate time elapsed since time to, t — to, and the elapsed time as computed by 
counting oscillator cycles, t* — to, is given by {t — to) — {t* — to) = t — t* = At* , where At* is the oscillator clock bias 
(error) that has accrued since the initial time to- So the time kept by the oscillator, t* , and the coordinate time, t, 
are related by 

t* ^t- At* (6) 

In the simplified model of carrier phase measurement Q , the carrier phase measured by a receiver is taken to be 
the phase difference between the incoming wave (at the antenna phase center) and that of the local oscillator, given 
by 

A0y (t) = (j),j{t) - (j)osc{t) = k\r,{t)- Rj(t - t)| - w At* + 27rmy (7) 

where the integer rriij is needed since the actual phase measurement can only be done up to an integer number 
of cycles. For example, if the geometric range between receiver and satellite is increased by an integer number of 
wavelengths, the measured phase difference, as given in Eq. (^, would be the same due to a different integer. 

By taking the origin of coordinates near the user platform, and taking differences of the measured quantities in 
Eq. (0), the clock bias At* can be eliminated from the unknowns, leaving only geometric quantities, such as the baseline 
vectors between antennas, r; — Vk, a unit vector to the satellite, Rj/|Rj|, and unknown integers. For example, for 
different antennae i, k tracking the same satellite j, the difference of phases given by Eq. (^ is 

A0y (t) - A(l)kj{t) = k (In - Rjl - \rk - Rjl) + 27r(my - nikj) (8) 
« fc |(rfe ~ r^) • + 27r(my - mfcj) (9) 

where |ri|, |rfc| << |Rj|. In the simplest case, to determine the platform attitude, we need three antennae to track 
two satellites. Two unknown baseline vectors, say (ri — r2) and (r2 — rs), uniquely determine the plane (attitude) of 
the platform. Each of these baseline vectors has three unknown components, so there are six unknowns. Therefore, 
tracking two satellites, j = 1, 2, with each of the three antennas, i = f , 2, 3, results in six measurements A0y(i). The 
six difference equations given by Eq. (^ can then be solved (up to integers) for the six vector components of the two 
baseline vectors, — r^, for say i — 1 and fc = 1, 2. 

The standard technique 01 for attitude determination, described in Eq. (|l|)-(^, uses a measurement model of 
carrier phase given by Eq. (Q), which essentially treats the electromagnetic field as a scalar field, and neglects its 
vector nature. The electromagnetic vector field carries additional orientation information that can be exploited to 
determine attitude of a user platform. Recently, it was shown that, using a short-dipole receiving antenna to track 
signals from six GPS satellites, two angles of attitude can be determined i.e., the orientation of the antenna can 
be determined up to rotations about its axis. The method used to derive these results was limited to a short-dipole 
type of receiving antenna, and to a determination of only two of the three Euler angles of attitude ||,|^ . 

In what follows, I show that using a single receiving antenna of arbitrary type, and tracking the carrier phase 
from seven satellites, the complete 3-axis attitude (in addition to position and time) can be determined. The results 
here are based on a model of the open-circuit voltage induced in the receiving antenna. The model explicitly takes 
into account the vector nature of the electromagnetic field that is broadcast by the satellite. The model shows that 
the phase of the voltage induced in the receiving antenna contains an additional phase-shift Ss{t), see Eq. (^2|), that 
depends on the relative orientation of the transmitter and receiver antennas. Consequently, given the position and 
orientation of the transmitting satellite antennae, the orientation (3-axis attitude) of the receiving antennae can be 
computed. 
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III. OPEN-CIRCUIT ANTENNA VOLTAGE 



The basis of the single-antenna method of attitude determination is the open-circuit voltage Vs{t) induced in a 
receiving antenna by the radiation field from satellite s: 

y,(t) = h(n,)-E,(R„,t) (10) 

where Es(Ro,t) is the electric radiation field at the receiver antenna located at position Rq at reception time t and 
h(ns) is the receiving antenna vector effective length (or effective height), which characterizes the antenna. The unit 
vector Us specifies the direction of field propagation, as seen in the comoving frame of the receiving antenna. The 
origin of this comoving system of coordinates is centered at the receiving antenna, which is at point O. See Fig. 1. In 
Eq. (p^), both the electric field, 'Es{'Ro,t), and the vector effective height, h(ns), are vectors in the comoving frame 
of reference of the receiver. The vector is given in terms of the polar coordinate angles 9s and 0s, in the receiving 
antenna's comoving frame: 

rig = — [sin 9s cos <j)s ai + sin 9s sin <j)s sl2 + cos 9s sls] (11) 

Here, ai, a2, as are the orthonormal unit vectors in the comoving system of coordinates that is centered at the 
receiving antenna. Neglecting relativistic aberration, this direction of propagation can be specified by a unit vector 
Us, pointing from the satellite to the receiver^ in the ECI coordinate system |l7|l . 

The vector effective height, h(ns) in Eq. dlQ), is a two-dimensional complex vector field on the unit sphere sur- 
rounding the receiving antenna that describes the receiving properties of an arbitrary antenna. In general, the field 
ii{9, (j)) is specified on the unit sphere by two complex functions, hg{9, (p) and h^{9, (p), 

hin) = hei9,(l))ag + h49,^)a^ (12) 

where 9 and (f) are the polar angles (coordinates on the unit sphere) , and a.g and a^ are orthonormal basis vectors on 
this sphere. The origin of coordinates is co-located with the receiving antenna. In practice, the functions hg(9, (f) and 
hif,(9, (j)) must be found by computation, or, by experiment. In the case of experiment, these two functions are found 
by mapping the far-field pattern of the receiving antenna, when it is used in the transmission mode. Then, using the 
equation for the far-field 11 



= Ee{9, 0) ae -f E^{9, cjy) a^ = -j^^ h(0, 0) e--'"('='^-*) (13) 

where is the current through the antenna, e is the permittivity, fi is the permeability of the surrounding medium, 
and r is the distance from the antenna, the components of the vector effective height, hg{9^(j)) and h^{9^4i), can 
be determined. Note that this procedure is well-defined for any type of receiving antenna. For a discussion of 
meas urements an d theory of elliptically polarized fields, see Ref. |^,^ 10 1^ and the coordinated series of papers in 



Rcf. |Jjl9[|20|j21|l. 

The electromagnetic field broadcast by the satellite, at the position of the receiver, is a far-field. In this far-field 
region, the general form of the electromagnetic field at receiver position Ro at time t is given by the real part of 

E(Ro,i) - [u£;„(e,$)e^"(®'*) -f v£;,e^-(®'*)]^^e^"* (14) 

r 

where the vectors u and v are real orthogonal unit vectors that specify the polarization of the field and the range, 
r = |Ro — Rt|, is the distance between reception event at the receiver (i,Ro) and emmission event at the satellite 
(tTjRr). The quantities Eu and Ey are real amplitudes that depend on polar angles Q and $, measured from the 
origin at point T at the satellite antenna. See Fig. 1. Each polarization can have an arbitrary phase and ipy. The 
field given by Eq. (|l^) is a general elliptically polarized electric field. 

As an example of a specific case of Eq. (|lj), consider a GPS satellite, where the field near the central beam axis is 
a right circularly polarized field, given by 



e 



-jkr 



E(Ro,t) = [u-jv]S(e,$) e^"* (15) 

r 

where E{Q, $) is a real amplitude. If the field broadcast by the satellite is not a right circularly polarized field, then 
Eq. (|l^) still applies, however, the vectors u and v will have complex phases, u = UoC-''^" and u = Voe^"^'' , where Uq 
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and Vo are real orthogonal unit vectors. Consequently, no generalization is lost by using Eq. ( |15D for the far-field of 
the satellite. This representation is used below. 

Strictly, E(Ro,<) in Eq. ( p^ is the electric field in the comoving frame of the satellite, while Es(Ro,i) in Eq. ( |lO| ) 
is the field in the comoving frame of the receiver. The difference between these two fields is a quantity of first order 
in u/c, where v is the speed of the satellite in the frame of reference of the receiver. In what follows, I drop these 
first order terms in v/c, taking Es(Ro,t) w E(Ro,<). Furthermore, the vectors u and v in Eq. ( |l5|) are taken at the 
retarded time t — rs/c, which is the time at which the signal left the satellite. 



Using the far- field of the satellite given by Eq. (15), the open-circuit voltage induced in the receiver antenna due 



to the electromagnetic field from satellite s, given in Eq. (10), can be written as (summation over repeated indices): 



Vs{t) = h,{n,{t)) D,k{t) R^k{t~^)-jR2k{t~~-) Bins) e^'^' (16) 



c c 



e 



-jkr 



r 



Here, r\.g{t) is the unit vector from satellite s at position Rs(i — ^g/ c) at emmision event time t — r^/c, to the receiver 



at position Ro(t) at reception event time t, given by |17 



Ro(t) - RJt- ^) 

where the positions of the receiver Ro(i) and satellite Rs(^)i ^re specified by their vector components in an inertial 
Cartesian system of coordinates, such as ECI coordinates, by 

Ro(t)=<(t)e, (18) 
B.s{t)^xl{t)e, (19) 

and Gi, i = 1, 2, 3 are the Cartesian basis vectors in the ECI (or other) inertial coodinates, xl,{t) and xl{t), i — 1,2, 3, 
are the time-dependent coordinates of receiver and satellite, respectively, and summation is implied when repeated 
indices occur, unless stated otherwise. The range from the satellite s emmision event to the reception event at the 
receiver, in ECI inertial coordinates, is given implicitly by 

= |Ro(t)-R,(i--)| (20) 
c 

The quantities hi{ns{t)) appearing in Eq. ( |l6| ) are the components of the vector effective height of the receiving 
antenna projected on the ECI Cartesian coordinate basis vectors e,; and are given in terms of the spherical basis 
components in Eq. (|l2|), by 



h(n) = hg{e, (j)) ae + h^{9, 0) = hi{n) ei -I- /i2(n) e2 + /i3(n) 63 (21) 

where 

hi{n) = hg{9, (j)) cos cos — h^{9, cj)) sincj) (22) 

/i2(n) = hg{9, (j)) cosOsintj) + h^{9, (f>) cos0 (23) 

hsin) ^ ~h0{9,(l))sm9 (24) 

These components depend on the direction (9, (j)) from which the radiation from satellite s is incident on the receiving 
antenna. This direction can be specified by components of the unit vector ns{t), defined in Eq. ([l7|). 

In Eq. (p^, the matrix elements Dik{t) specify the attitude of the receiver platform at reception time t. The 
attitude is specified by giving the inner product of the unit vectors a^, i = 1,2, 3, of the receiver comoving coodinates 
with the ECI Cartesian coordinate basis vectors e^: 

D,k{t) = &^{t) ■ e^ (25) 

where the vectors ai(f) are time-dependent, since these vectors correspond to the moving body axes of the receiver. 
The receiver platform attitude matrix Dik(t) — Dik{a{t), P{t),j(t)) has three independent parameters. This attitude 
matrix can be conveniently represented in terms of three time-dependent Euler angles, a, (3, and 7, by |22[ 

cos (3 cos a cos 7 — sin a sin 7 cos (3 sin a cos 7 -I- cos a sin 7 — sin /3 cos 7 
Dik{t) — I — cos /3 cos a sin 7 — sin a cos 7 — cos /? sin a sin 7 -I- cos a cos 7 sin /3 sin 7 | (26) 
sin /3 cos a sin [i sin a cos fi 
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The attitude matrix Dik{t) completely specifies the orientation of the basis vectors of the receiver platform, a^, in 
terms of the ECI basis vectors e^, see Figure 1. 

In Eq. ([l6|), the quantities Rik{t — rs/c) and R2k(i ~ ^s/ c.) specify the satellite antenna polarization vectors, u and 
V, at the time of transmission, t — Vsjc^ with respect to the ECI Cartesian coordinate basis vectors e^: 



T 

\i(t ^) = (u • e„i) e„, = R\„i(t - 

c 


C 


(27) 


v(t ^) = (v • e„) e,„ = R2rn{t - 

c 


C 


(28) 
(29) 


In Eq. ([l^), the remaining quantities to be defined are 






c c 




(30) 



where = N^i + jNs2 is the effective index of refraction (real and imaginary) of the medium through which the 
electromagnetic wave travels, from satellite transmitter to the user's receiver. The frequency / is the L-band (or 
other) carrier frequency. 

The voltage in the receiving antenna due to signal from satellite s, given in Eq. (Iq), can be written as 



V,it) = V, e-^'^ E{ns) ^^^e^'"^* (31) 
where Vg is a real amplitude and Sg is a real phase angle, specified by (sums on repeated indices i and k): 



Ssit) = -ATg^h^{nsit))D^kit) 



Here, for real numbers a and (f), the function Arg(ae-''^) = (p mod 2tt. 

The actual voltage in the receiving antenna is the real part of Eq. (Hq) , which can be written as 



(32) 



Re V; (t) = vs cos ips (t) (33) 

where 

ipsit) = k,r,it)-ujt + S,{t) (34) 

The phase of the voltage in the receiving antenna, '0s(^) given in Eq. (Q), should be compared with the phase of 
the scalar field, given in Eq. (||). The vector nature of the electromagnetic field leads to an additional phase shift 
6s{t) in Eq. (34). This phase shift is a function of the relative orientation of the satellite and receiver antenna. The 
single- antenna method of determining attitude exploits the dependence of the phase 5s{t) on the relative orientation 
of the satellite and receiver antennas. Note that Ss{t) in Eq.(^) depends on receiver antenna parameters hi{ns), 
which are defined in Eq. (p^)-(p^). 

The model for the phase of the induced voltage in a receiving antenna, given in Eq. (^^ , is more accurate than the 
standard model in Eq. (^, because it takes into account the orientation of the receiving antenna in the electromagnetic 
field. 



IV. IMPROVED CARRIER PHASE MEASUREMENT MODEL 

The satellite receiver has an electronic oscillator, whose phase is given by Eq. ([4|). The measured quantity is the 
phase difference between the local oscillator and the phase of the voltage in Eq. (pif): 

V-s (t) - V'osc(i) = ks rs (t) + ujAt* + 6s {t) + 2TTm (35) 

where m is an integer, since the relative phase can be measured only up to integer multiples of 27r. For comparison 
purposes, Eq. ( |35| ) can be put into a notation similar to Eq. (^: 

tpisit) - i>osc{t) = hs ris{t) + LoAt* + 5is{t) + 2TTmis (36) 

where the subscript i has been added to index the receiving antenna, and the subscript s still refers to transmitting 
satellite s. Comparison of Eq. (^) with Eq. (^) and the previous definitions for 5s{t) makes the meaning of the 
terms in Eq. ( |3^ ) clear. 
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V. ATTITUDE DETERMINATION AND NAVIGATION 



A. Multiple- Antenna Attitude Determination 

For the case of an attitude determination system that uses multiple antennas, Eq. ( p6[ ) can now be used to form 
differences between two antennas and one satellite, such as is done in standard carrier phase tracking (see Eq.(||)): 



A'0is(i) - AV'fcs(i) = {hsTisit) - kksrksit)) + 27r(mis - mks) + Sis{t) - dks{t) (37) 



^ { - • 1^ + 27r(m,, - mus) + 5is{t) - (38) 



where are the positions of the receiving antennas at reception time t. The measured phase, as given by Eq. ( |3j 
contains an additional difference of phase shifts, 5is{t) — Sks{t), which depends on two attitude matrices, Dis and 
Dfcs, one for each of the two antennas mounted on the user platform. Compare Eq.(38) with Eq.(||). When the two 



antennas are identical, and they are mounted on a rigid platform, and they are aligned, then these two matrices are 
the same. The difference of phase shifts Sis{t) — Sks{t) is essentially zero. Alternatively, if these conditions are not 
satisfied, then Eq.(^ provides the generahzation of Eq.(^), and includes the effects of flexing of the platform and 
misalignment of the antennas. 



B. Single- Antenna Attitude Determination 



The orientation dependence of the phase shift 5i [t] in Eq. ( p2[ ) can be exploited to permit determination of attitude 
(by determining a, (3, and 7) from carrier phase measurements made by a single antenna. Consider a receiver whose 
antenna has a known orientationa and position at some initial time to- The phase change of the oscillator in the 
receiver, since the initial time to can be represented by 

AV'o.c(0 = ~^{t* ~ to) (39) 

where t is the GPS time (or coordinate time in the ECI (or other) frame of reference) that is associated with a time 
t* on the clock that counts oscillator cycles ||l^. The oscillator in the user's receiver is assumed to be imperfect. The 
imperfection of the oscillator can be represented by the fact that the oscillator does not keep time according to ECI 
coordinate time. The oscillator will then have a running phase given by 

4>osc{t) = -tot* (40) 

where the time scale t* is not the same as the ECI coordinate time t. I assume that at the initial time to the oscillator 
has the correct time t* = to. If we counted the cycles of the oscillator since the initial time to, and multiplied by the 
period of the oscillator, 27t/uj, then the elapsed time, according to the oscillator is t* — to- The difference between the 
real ECI coordinate time elapsed, t ^ to, and the elapsed time as computed by counting oscillator cycles, t* — to, is 
given by {t — to) — {t* — to) — t — t* = At* , where At* is the oscillator clock bias that has accrued since the initial 
time. So the time kept by the oscillator, t* and the coordinate time, t, are related by 

t* =t- At* (41) 

Using Eq. (|l]) in Eq. §) gives the oscillator phase change since the initial time to in terms of ECI coordinate time 
and the oscillator clock bias 



A^Poscit) = -ujit - to ~ At*) (42) 

The phase of the voltage (due to satellite s) in the receiver antenna is given by Eq. (]3^) . Since the initial time to, 
the phase of the voltage in the antenna changed by 

At^sit) = A{t) - Mto) = h (rsit) - rs{to)) - uj{t - to) + 5s{t) - 5,{to) (43) 

By taking the difference of the phase changes that have accrued on the receiver's oscillator (see Eq. (^2|)) and the 
carrier phase from satellite s (see Eq. (^)), I get the single antenna carrier phase-change equation 

A^',(t) = ATPs{t) - A^oscit) = ks (r,(<) - Tsito)) -LuAt* + Ssit) - Ssito) (44) 
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Equation (^4|) gives the difference of the accrued phase changes in the oscillator and in the voltage of the receiver 
(due to signal of satellite s). This is the quantity that is measured by satellite tracking receivers and is known as a 
measurement of carrier phase or accumulated delta-range. Equation (Q), together with the above definitions, defines 
an algorithm for single-antenna determination of position, time, and attitude by satellite navigation. 

The measured phase change given by Eq. (|4^) depends on seven unknown parameters and several assumed-known 
parameters: 

A^,{x'^,At*,a, (3, r, x'l.Rik, R2k,ao, Polo) (45) 

The unknown parameters at time t are the receiver antenna position, x^, fc = 1, 2, 3, receiver oscillator phase correction 
LolSi* (or equivalently, the receiver clock bias Ai*), and the three Euler angles, a, j3, 7, which define the attitude 
matrix of the receiver antenna. The parameters that are assumed known functions of time are the satellite position 
coordinates or satellite ephemeris, x^', k = 1,2,3, and the satellite attitude parameters Rik, i?2fc, fc = 1,2,3. I will 
collectively refer to a;^, Rik and i?2fe, as the satellite enhanced ephemeris. In addition, the receiver attitude at the 
initial time to must be known through, for example, the Euler angles ao, /3o, and 70. 

There are several modes in which a user may exploit the general single antenna attitude determination algorithm 
described above. The choice of mode depends on other systems that are available to the user, and what information 
the user requires. Perhaps the most basic assumption is that the user has no additional systems for navigation. In 
this case, the user receiver determines the needed parameters at the initial time as well as the satellite enhanced 
ephemeris, and then tracks the carrier phase continuously of seven satellites. At a given time t, the seven equations 
in Eq. (Q), for s = 1, • • • , 7, are a closed system of equations that can be solved for the seven parameters of receiver 
position, Xq, fc = 1, 2, 3, receiver clock bias At*, and the complete three-axis attitude given by the three Euler angles 
a, /3, and 7. The requirement that the user's receiver track seven satellites is not as severe as it may seem, because, 
for example, in the case of GPS, there are plans to provide additional satellites to the number that are currently 
operational. 

Another mode of navigation is to carry an oscillator that is sufficiently stable over the intended time interval of 
navigation. For example, in the case of a UAV, the time of flight may be 10 minutes. An oscillator accurate to one 
part in 10^ would maintain a phase within 

LoAt* = 27r(1200.s"i)((600.s)(10~^)) = 0.0045 radian « 0.26° (46) 

If the oscillator is accurate enough, then the parameter A*< does not have to be determined, and it is sufficient to 
track six satellites, to compute the three position coordinates and the three Euler angles of attitude. 

Still another mode of navigation is to provide the user with compass heading information. The ability to determine 
three axis attitude allows the extraction of the compass heading associated with a particular axis of the receiver 
platform. This requires a trivial transformation from ECI coordinates and time to topocentric coordinates. 



VI. SUMMARY 



Current technology of attitude determination by satellite navigation requires the use of multiple, well-separated 
antennas. In this paper, I have described the theory and the algorithm for determining the position, time and three 
angles of attitude (antenna orientation), when tracking the carrier phase of seven satellites using a single antenna. 
The method is based on the phase change induced in the open circuit voltage in an antenna due to the signal of a 
satellite that is transmitting from a known location and having a known orientation. The dependence of the phase 
shift on receiving antenna orientation is given in Eq. (|3^). 

The attitude determination algorithm described here can be applied by current GPS receivers with essentially a 
software addition, since these receivers routinely track the carrier phase (or delta range) of the satellites. Furthermore, 
the orientation of the GPS antennas is known quite accurately because the satellite main beam is pointed toward the 
center of the Earth . 

If a higher accuracy of attitude is required, the satellite broadcast ephemeris can be enhanced to include precise 
information on the satellite antenna orientations. This enhanced ephemeris could be implemented in a future version 
of the GPS, or some other satellite system. Such an architecture would make the satellite system capable of providing 
a user with time, position, and orientation. 

For user applications, a single-antenna (rather than multiple-antenna ) attitude determination system is preferable 
because it allows miniaturization of the receiver system and is potentially less expensive to manufacture and to 
install. An example where single- antenna attitude determination is preferable is in micro air vehicles (MAV), nano 
satellites |l^ and in general aviation aircraft. Currently, the U.S. fleet of aircraft is transitioning from ground- 
based navigation methods to satellite-based (mostly GPS) navigation. Attitude information is needed by a pilot 
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whenever aircraft are flying by reference to instruments, such as in instrument meteorological conditions (IMC). 
Many general-aviation aircraft already have single-antenna GPS receivers. Therefore, it may be possible to extract 
attitude information, using the single-antenna attitude determination method presented here, from GPS carrier phase 
tracking data in already-installed GPS receivers, with only a software change and addition of an attitude output 
channel. 

Information on position, time and attitude of a platform is equivalent to knowing the compass heading. Conse- 
quently, another application of attitude determination is to replace the wet compass. Satellite navigation receivers 
are being miniaturized and put into a wrist watch ||l5[| . The potential for miniaturizing the user receiver by using 
the single-antenna attitude determination method may lead to production of a satellite-based compass that can be 
carried by a person in a watch-sized instrument. 

Finally, the algorithm presented here can be applied to GPS signal simulators, to simulate a real signal seen by a 
satellite receiver that is rotating as well as translating. 
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FIG. 1. The Earth-centered coordinates with basis vectors e^, i = 1, 2,3, is shown together with the comoving system of 
coordinates y' centered at the receiver (point O) and the satellite comoving system of coordinates centered at point T. The 
unit vector n points from the emission event at the satellite at time t — r^/c to the reception event at the receiver at time t. 
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